Comparing mesoscale chemistry-transport model and 
remote-sensed Aerosol Optical Depth 

^ ^ C. Carnevale^, G. Mannarini '^'^'*, E. Pisoni^, M. Volta'^ 

O ' 

fvq . "• Dipartimento di Ingegneria deU'Informazione, Universita degli Studi di Brescia via 

^^ ■ Branze 38, 1-25123 ITALY 

^" 

a^ 

Abstract 

p ,■ A comparison of modeled and observed Aerosol Optical Depth (AOD) is 



O 






presented. 3D Eulerian multiphase chemistry-transport model TCAM is 



c^ , employed for simulating AOD at mesoscale. MODIS satellite sensor and 

Y^ \ AERONET photometer AOD are used for comparing spatial patterns and 

•^ I temporal timeseries. 
^' TCAM simulations for year 2004 over a domain containing Po- Valley and 



nearly whole Northern Italy are employed. 

For the computation of AOD, a configuration of external mixing of the 
-^ , chemical species is considered. Furthermore, a parametrization of the effect 

>■ I of moisture affecting both aerosol size and composition is used. 

S:: ' An analysis of the contributions of the granulometric classes to the ex- 

00 ■ tinction coefficient reveals the dominant role of the inorganic compounds of 

submicron size. For the analysis of spatial patterns, summer and winter 
^ \ case study are considered. TCAM AOD reproduces spatial patterns similar 

Q \ to those retrieved from space, but AOD values are generally smaller by an 

order of magnitude. However, accounting also for moisture, TCAM AOD 
significantly increases to values of the same magnitude of the observed ones. 
^ . The temporal performance of model AOD is tested in correspondence of 

H \ AERONET site "Venise" and some temporal structures are reproduced. 

The results suggest encouraging perspectives in view of satellite AOD 
assimilation also at the mesoscale and not only at the global and regional 
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scale as already tested in the literature. 
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1. Introduction 

Chemistry- Transport Models (CTMs) often are the sole valuable tool for 
investigating the impact of pollutant emissions and air quality strategies at 
the mesoscale. This is due to the nonlinear nature of the process o f poUu - 
tant chemical transformation, diffusion and advection ISchaap et al.l ( 120041 ) . 
CTMs provide both gas and aerosol phase state variables, which either ex- 
actly represent the same chemical quantity measured at the monitoring sites 
(such as O3 or NO^;) or can be aggregated for computing monitored quanti- 
ties (such as PMIO). While CTM performance in reproducin g observed gas 



phase quantities is fair, PMIO simulation skill is generally poor lVautard et al 
(120071 ) ■ Among the reasons for CTM weak p erformance, i t is p ossible to 



list: uncertainties in the emission inventories iHenze et al.l ( l2009l ). not de- 
tailed representatio n of aerosol phase ch emistry, in particular the formation 
of secondary aerosol IStanier et al.l (120081 ). as well as difficulties in comparing 
point-like experimental information from mo nitoring sites to spatially aver- 
aged values from models lVautard et al.l (|2007| ). As long as these issues remain 
unsolved, CTM aerosol simulation and prediction power relies on data assim- 
ilation, which makes use of auxiliary observations (e.g.: air quality stations, 
remote sensing data) in order to reduce the uncertainties in input d ata, such 



as initial conditions or boundary conditions iTombette et al.l ( 120091 ) 



Spatially and temporally homogeneous observations on the mesoscale can 
be provided by means of optical instruments for remote sensing of aerosols. 
Aerosol Optical Depth (AOD), which is the column integrated extinction 
coefficient, represents the optical effect of the total aerosol load in the atmo- 
sphere. There are both ground based stations as well as airborne measure- 
ments of AOD. 

A n AOD output has bee r i computed from vario us CTM s, both at the 
globa 



Morcrette et al.l(l2009l) : lBenedetti et al.l (20091 ). regional lKoelemeiier et al. 



(|2006h : lAdhikarv et al.l(l2007l):IChoi et al.l (l2009l) . an d mesoscale JHodzic et al. 



(l2004j ): lde Meii et al.l (l2007f ): IChapman et al.l (120091). This brief survey is de- 
liberately restricted to models run at the mesoscale. 



CHIMERE model was run over the Paris greater metropolitan area with 
a horizontal spatial r esolution of 6 km an d 20 vertical layers up to the pres- 
sure level of 500 hPa iHodzic et al.l ( 12004 ) . This model employes a sectional 
approach with 6 bins between 10 nm and 40 //m. The model delivers both 
a dry and a wet particle radius, i.e. a radius accounting for the hygroscopic 
growth. The wet radius is obtained from a parametrization which depends 
on relative humidity (RH) only. Vertical profiles of the aerosol backscat- 
tering coefficient were computed and compared to Lidar measurements. The 
average complex refractive index and the Single Scattering Albedo were com- 
pared to Sun photometer measurements. It was found that the aerosol load 
is typically underestimated by the model and that a prevalence of submicron 
sized particles is simulated. 

In another mesoscale modeling exercise, AO D was comput e d by T APOM 



model over a portion of the Po- Valley (Italy) Ide Meij et al.l ( 120071 ). Hori- 



zontal spatial resolution down to 5 km and 21 vertical layers up to 6000 m 
a.s.l were employed. The model considers a sectional approach with 4 fixed 
size bins ranging from to 10 /im. The activity coefficients of atmospheric 
aerosols are assumed to be in equilibrium with the environment, governed by 
RH. Model AOD was compared to both MODIS and MISR satellite obser- 
vations as well as to AERONET and CHARM Sun photometers data. The 
model was found to underestimate observations. 

Also the WRF-chem mo del has bee n deve l oped for comput i ng ae rosol 
chemistry and aerosol optics iFast et al.l ( 120051 ): IChapman et al.l (120091 ). A 
domain comprising western Pennsylvania was modeled with 2 km resolution. 
Vertically, the domain extended 57 nodes in the vertical, from the surface 
to 100 hPa, with finer resolution near the surface. The model was run with 



8 dirn ensional sections, ranging from 20 nm to 5 fim radius IChapman et al 
(120091 ). Thanks to a specific module for calculating the activity coefficients 
of various electrolytes in multicomponent aqueous atmospheric aerosols and 
a solid-liquid equilibrium solver, the model is able to describe aerosol water 
uptake. However, particle transfer among sections as a consequence of water 
uptake is forbidden. Model AOD was compared to radiometer observations 
capturing the diurnal cycle in sequence of a few days, although predicted 
AOD on a day after precipitation events is clearly too high. 

In the present paper, an optical module is presented which works on top 
of the aerosol simulations from TCAM (Transport and Chemical Aerosol 
Model). TCAM is a multiphase three-dimensional Eulerian model. A pecu- 
liarity of the model is that aerosol particles are described by 10 granulometric 



classes dynamically set by the model, allowing a detailed desc ription of the 
chemical reactions affecting particle size lCarnevale et al.l (l2006l ). The domain 
of investigation is Po- Valley, in Italy. Northern Italy is a densely inhabited 
and industrialized area, where high anthropogenic emissions and frequent 
stagnating meteorological conditions regularly cause high PMIO levels. For 
the computation of AOD, an external mixing approach of the aerosol species 
has been chosen. Hygroscopicity is accounted for in a parametric way. For 
the sake of making an homogeneous comparison with satellite AOD prod- 
ucts, a cloud mask is applied to the model outputs. At the AERONET site 
"Venise", both model and satellite AOD are compared to the reference value 
from the photometers. 

The paper contains a detailed description of the optical module for the 
computation of AOD from TCAM aerosol fields Sect. ([2]) and a case study, 
including comparisons to observations, which is presented in Sect. ([3]). A 
combined Discussion and Conclusions section is also provided Sect. (jl]). 



2. TCAM 

2.1. Aerosol Module 

TCAM is a part of t he Gas Aerosol Modeling Evaluation System (GAMES) 
Volt a and Finzil (120061 ) which also includes the mete orological pre-processo r 
PROMETEO and the emission processor POEM-PM lCarnevale et al.l fl2006h . 

Within TCAM the chemical and physical dynamics of aerosol particles 
of different sizes is modeled. The aerosol module considers 10 size bins (usu- 
ally ranging between lOnm and 30/xm), whose bounds can vary during the 
simulation, according to the gas-to-particle and growth phenomena, and ac- 
cording to the concentration of each compound. This detailed description 
of the granulometric classes within TCAM is an important difference with 
respect to other CTMs. 

The aerosol module implemented in TCAM is coupled with the CO- 
COH97 gas phase chemical mechanism. TCAM describes the most relevant 
aerosol processes: the condensation, the evaporation, the nucleation of H2SO4 
and the aqueous oxidation of SO2. 

Water is assumed to be always in equilibrium between the gas and the 
aerosol phases, i.e., water activity of an aqueous aerosol is equal to the am- 
bient RH. Thus, water is not a transported quantity: it just mirrors the 
met eorological input. This assumptiori is als o made by other model develop- 



ers f ide Meii et al.l . 120071 IWexler et al.l . 119941 ) 



More technical detail s abou t the model and its performances may be 



found in ICarnevale et al.l (120081 ) . 



2.2. Optical module 

The optical module takes as an input the aerosol mass concentrations 
resulting from TCAM simulations and delivers the associated AOD. The 
approach is based on the computation of the extinction coefficient a, from 
which AOD is obtained after integration over the optical path, 

rhii 

AOD= / a{z)dz (1) 

Jhi 

Aerosol particles are considered as externally mixed. Thereafter, a results 
from linear superposition of the contributions a^*-* of different aerosol species. 
Furthermore, each of the hk model size classes has to be added up for ob- 
taining a^^\ This leads to the double summation: 



ns nji 

a = 

s=l k=l 



EE4" (2) 



Individual addends in the extinction coefficie nt are given by the numbe r 



concentration times the extinction cross section iBohren and Huffman! ( 119831 ). 

ai'^ = Nt^ aiiin) (3) 

(s) 

The extinction cross section a^J^ is the geometrical cross section of the par- 
ticles (assumed to be spherical) multiplied by the extinction efficiency. The 
model concentrations are considered as related to homogeneous spheres of 
radius r^. Thus, Mie's extinction efficiency Qext comes into play: 

aiiin) = irrl Q,,,{n; A = 550nm; m^'^) (4) 

Here the dependence on the complex refractive index m*^*^ = n^'^^ — in^^'' of 
specie s is highlighted and the in-air light wavelength A is kept fixed. The 
extinction efficiency plays a crucial role in the calculation of AOD since, for 
Re{m) < 2, its value varies from exactly zero to about 6. 



2. 3. Effect of moisture 

Some of the atmospheric aerosol species are hydrophihc: above a certain 
RH they are no more in the sohd state but rather form an aqueous solu- 
tion. Depending on chemical composition, there is eithe r a continuous or an 
abrup t (deliquescence) transition to the aqueous solution ISeinfeld and Pandis 
( 1l998l ). As a result of water uptake, both particle density and refractive in- 
dex are driven towards liquid water density and refractive index respectively. 
This affects particle optical properties twofold. On the one hand, the aerosol 
size distribution is shifted towards larger radii as a result of hygroscopic 
growth. On the other hand, the extinction efficiency curve as a function of 
radius Qext{r) is modified. The position of its maximum is shifted towards 
a larger value because wet refractive index is generally lower than dry one. 
Which of either effects prevails, extinction enhancing or quenching, depends 
on the actual size of peak shift and hygroscopic growth. 

TCAM dry particle radii r^o are rescaled by a RH-dependent factor. 



rk = Tko 



■P^'\RH) 



(5) 



Such wet radius r^ is then employed in Eq. (jl]) for the computation of the 
extinction cross section. At this place it is worth noting that the aerosol 
number concentrations A*"^** are obtained from dry radius Vko'- 



N, 



^'^ - MtU P^'^^nrl, 



(6) 



with the mass densities p^^-* as in Tab.[T] NJf are defined on a 3-dimensional 
grid given by the cartesian product of a 2-dimensional horizontal grid with 
njy=ll vertical layers of TCAM. 



I3'^^\RH) is taken to be a fit of Tang's calculations iTangI ( 1l996l ). which 
are based on the single particle levitation technique measurements. In doing 
so, any hysteresis in the growth-evaporation cycle of hydrophihc particles 
is neglected. A sin g le growth curved has be en used also by other authors 
Hodzic et all f l2004h : [Koelemeiier et"aD fl2006[ ). 

As a second step for accounting for moisture, wet complex refractive index 
m = n — in is computed from a weighted average of dry aerosol index rriQ and 



liquid water index m^, weights being the volume fractions in the mixture 

nioPo + m„P„ 



m{RH) 
p 



Po + Pw 



(7) 



1 



[^^'>{RH)Y 



1 



(The weights Po,w are expressed in units of dry particle volume). Volume 
weighting is chosen, being the simplest of three weighting procedures which 
are nearly equivale nt when applied to the computation of refractive index 
Lesins et all f l2002h . 



A fair fit of the growth factor displayed in iTangl (jl996[ ) for a group of 
inorganic compounds is the following function: 



^'^'\RH) = 1 + ^ 



is) 



-RH 



-iRH) 



12 



(8) 



Eq. dH]) is used for RH< 100%. In case of supersaturation, capping of hy- 
groscopic effects is applied, i.e. f^^^^RH > 100%) = 1 + 87'^''^. The sole 
adjustable parameter in Eq. ([H]) is 7*^**'*, whose values are listed in Tab.[Tl 

The hygroscopic properties of organic parti cles emitted d u ring l aboratory 
biomass burning experiments were studied in iPetters et al.l ( 120091 ). The re- 
sulting growth factor strongly depends on the investigated specie. However, 
TCAM organic aerosols include alkanes, aromatic compounds, olefins, and 
cresol. None of them is hydrophilic, since they are not able to bind water 
or (case of cresol) since they are not well soluble. Thus, 7^*'' = for such 
species. 

2.4- Cloud screening 

MODI S retrieval a r e not allowed over cloudy or foggy pixels by the AOD 
algorithm iLevy et al.l ( 120071 ). Conversely, model AOD is available even in 
case of total cloud coverage. Thus, for making satellite and model AOD 
more directly comparable, a cloud mask is applied to model AOD too. The 
diagnostic variable Cloud Water Content (CWC) is employed, which is ob- 
tained from MM5 outputs. The screening protocol consists in requiring null 
CWC over the whole vertical profile. The protocol ensures that all cloud 
contaminated pixels are screened out from the model. 



compound p 


g/cm^ 


n 


K 


Ref. 


7 


EC 


1.80 


1.95 


.790 


Bond and Berestrom f2006) 





OC 


1.50 


1.53 


.006 


Tsiearidis et al. f2008l 





sulfates 


1.60 


1.53 


.000 


de Meii et al. ('2007) 


1 


nitrates 


1.60 


1.53 


.000 


de Meii et al. f2007l 


.75 


ammonium 


1.60 


1.53 


.000 


de Meii et al. (2007) 


.75 


dust 


2.65 


1.50 


.033 


Reddy et al. (2005) 






Table 1: Mass density p, real n and imaginary k part of dry refractive index (at wavelength 
X=550 nm), hygroscopicity parameter 7 (see Eq. (E^), for each of the aerosol compounds 
considered in this work. 



3. The case study 



The method outlined in Sect. (12. 2 p has been applied to TCAM aerosol 
simulations from January 1 to December 31, 2004. The study area covers 
nearly whole Northern Italy, which is divided into 41x64 cells (10x10 km^ 
resolution) In the following this area is also termed " QUITSAT domain" . 

Simulations have been performed getting initial and boundary conditions 
by a nesting procedure fro m the results of CHIMERE simulations over Eu- 
rope |Sclunidt_et_alJ (1200 ll ). The emission fields were estimated by POEM- 
PM pre-processor starting from the regional emission inventory collected by 
the ISPRA (Istituto Superiore per la Protezione e la Ricerca Ambientale). 
The point sources and area emissions have been chemically and temporally 
adjusted to TCAM requirements and have been derived by CTN-ACE Ital- 
ian modeling intercomparison project . Meteorological fields have been cre- 
ated with MM5 model and adapted to TCAM through PROMETEO pre- 
processor. 

Hourly aerosol concentrations from TCAM QUITSAT v. 3 simulations 
have been employed as an input to the AOD computation. TCAM delivers 
the concentrations of 21 species, which have been aggregated into ns = Q 
chemical compounds. One of these is unbound or elemental carbon, "EC". 
Primary and secondary organic carbon species have been lumped into the 
"OC" class. Then, some inorganic electrolytes have been neglected for the 
sake of lighten AOD computation, since they exhibit very low concentrations 
in TCAM QUITS AT simulations: Th ey include sea-salt ionic components 
Na"^ and Cl~ (cp. Ide Meij et al.l (l2007l )). aqueous phase S02(aq), Il202(aq), 
03(aq), and protons H"*". Liquid water has not been considered as an aerosol 
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specie for its own. However, the role of water is accounted for through the 
parametrization of hygroscopic growth according to Eq. (jS]). Finally, one is 
left with the 6 aerosol compounds, as listed in Tab.[T] 

3.1. Size classes 

Each of the 6 chemical compounds is organized into ?7.i^=10 dimensional 
dimensional classes, whose size is different for each simulation pixel and for 
each time step, due to the specific fixed-moving algorithm for TCAM size 
distribution. In Fig.[T] domain mean size distributions averaged over time 
spans of one month are displayed. In these plots, the 6 chemical compounds 
are further lumped into an inorganic group (ammonium, nitrates, sulfates) 
and the remaining compounds (EC, OC, dust). Summer and winter typical 
situations are displayed. The 4 smallest dimensional classes (radii < 140nm) 
are always the most populated in terms of number concentration A^^. The size 
distribution of the inorganic group nearly follows a power-law distribution, 
while the distribution of the remaining compounds is peaked at a radius 
below 100 nm. This holds both in June and December, 2004. It is also 
seen that the smallest 3 classes are affected by the most significant changes 
between June and December, while the larger classes are rather insensitive to 
the choice of the month. The general trend is that larger radii are obtained 
in summer, independently of the chemical group. This behavior should be 
related to photo-induced aging of primary pollutants. 

Mie's extinction cross section aey±{r,\ = 550nm, m) with m = 1.53 — 
iO (1.95 — zO.79) for the inorganic (rest) group is also displayed in Fig.Hl It 
is obtained from a code by C Matzler (Universitat Bern), publicly avail- 
able on the web (www.iap.unibe.ch/publications/download/201/en/). The 
asymptotic behaviors of Mie's efficiency are exploited for lighten numerical 
computations. A look-up table of Qext{f, n, k) for logarithmically spaced radii 
r and imaginary parts k and for linearly spaced real parts of refractive index 
n is linearly interpolated at each model cell, at each size bin, at each time 
step. A look-up table appr oach to Qg ^ t is em ployed in the optical module of 



the WRF-chem model too lFast et al.l (12005f l. 



If the cross section ae^ti^k) gets multiplied by the aerosol number con- 
centration Nk, the class- resolved extinction coefficient ak is obtained. Since 
(7ext{fk) is (nearly) logarithmically increasing and Nk is logarithmically de- 
creasing for particle class-radius r^ > 80nm, the product of both has got a 
maximum value at intermediate values. Thus, the most significant contribu- 
tions to the total extinction coefficient result from particles whose radius is 
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Figure 1: TCAM number concentrations Nk (crosses); extinction cross section acxt{fk) 
and its geometrical limit (lines); extinction coefficient components ak (circles) for both dry 
(magenta) and wet (blue) case. Wet case is obtained considering RH=80%. Left panels: 
inorganic compounds (ammonium, nitrates, sulfates); Right panels: remaining compounds 
(EC, OC, dust). The remaining compounds are not hygroscopic, thus dry and wet curves 
overlay. For this group, the EC extinction cross section is displayed. Values in top (bottom) 
panels refer to the 3-dimensional domain average for the month June(December) 2004 
considering times between 9:00 and 1^:00 UTC 
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Figure 2: Temporal evolution between 9:00 and 1^:00 UTC of vertical profile of extinction 
coefficient a{t, z) for mean day in June 2004, accounting for moisture. Left(right) panel: 
model cell next to AERONET site Ispra(Venise) . 



below 0.5//in. Larger values of the extinction are obtained for the inorganic 
group, which is a consequence of larger population of this kind of aerosols 
according to TCAM. 

In Sect. (12. 3 p it has been mentioned that, in general, moisture may lead to 
both extinction enhancing or quenching. In Fig. [T] it is seen that in the actual 
case moisture results in enhanced extinction coefficient for every dimensional 
class. This is due to the fact that most of the aerosol population is submicron 
sized. 



3.2. Vertical structure 

The temporal evolution of the vertical resolved extinction coefficient a{z) 
is displayed in Fig. [2] for the mean day in June 2004. Data refer to the 
domain cells which are next to Ispra and Venise AERONET sites respectively. 
Largest a{z) values are obtained at heights z < 500m. For the highest 
layers, the extinction is about 1/10 of the maximum value assumed in each 
vertical profile. This means that the vertical modeling within TCAM is 
able to capture most of the optical extinction responsible for AOD (June 
usually is the month during which the insulation and thus the boundary 
layer height reach their maximum). In any case, a{z = hn) values for TCAM 
QUITSAT simulations are quite low and are comparable to the extinction 
resulting from molecular sc attering in a clean atmosphere 
Seinfeld and PandisI f tToish l 
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Figure 3: AOD maps for June 10, 2004 over the QUITSAT domain. Top-left panel: 
MODIS product "Corrected_OpticaLDepth_Land_And_Ocean" from both Terra and Aqua 
overpasses; Top-right panel: TCAM AOD averaged over the times corresponding to the 
satellite overpasses; Bottom-right: TCAM AOD accounting for moisture; Bottom-left: 
TCAM AOD accounting for moisture and for mask. Notice the vertical scale in MODIS 
AOD panels. 



3. 3. Comparisons of spatial patterns 

According to Eq. ([1]), TCAM AOD is obtained from integration of the ver- 
tical profiles of the extinction coefficient. Additionally, AOD can be obtained 
from satellite observations too. MODIS product " Corrected_OpticaLDepth_Land" 
is consi dered in th e follow ing, which is the quality reference for quantitative 
studies iLevy et al.l (120071 ) . This product is available at the wavelength of 550 
nm. Data from both Terra and Aqua satellite are considered. Daily mean 
AOD maps are obtained upon averaging over all scenes having finite spatial 
overlap with the QUITSAT domain. 

In the first panel of Fig.[3l such a map is shown for the day June 10, 2004. 
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This is a cloud-free day for most of the domain pixels, with exception of the 
alpine region and the western part of Tusco-Emilian Apennine. A strong 
aerosol signal arises from the eastern part of the domain, with AOD values 
close to 1 over Veneto and Friuli. Also, AOD ^ 0.4 is observed southern of the 
western Alps, which is nearly double as the value in the rest of western Po- 
Valley. Corresponding model AOD can be obtained averaging over all model 
maps whose simulation time corresponds within 30 minutes to the time of 
satellite overpasses. The results of this procedure are shown in the top-right 
panel of Fig.|3l In contrast to satelhte map, model map is defined for each 
domain cell. Its main feature is a hot spot extending from southern part of 
Veneto to eastern part of Emilia-Romagna. Then, an AOD arc southern of 
western Alps is present. AOD values are a factor of 2-5 lower than in the 
satellite map. However, the above described qualitative features of the map 
are present both in model and observations. 

The model output considered so far does not account for moisture. Its 
effect is computed as in Sect. (12. 3 p and can be appreciated in TCAM map in 
the bottom-right panel of Fig. [31 It is clearly seen that moisture leads in this 
case to AOD increase, which is now nearly double with respect to the dry 
map. The imp ortance of mo i sture for haze over Po- Valley has been already 



highlighted by lde Meij et al.l (120071 ). Furthermore, quite high AOD values on 
western Alps and Apennines cells are found. The origin of these high values 
should be related to Cloud Optical Depth. Indeed, in the bottom-left panel 
of Fig. El where cloud masking procedure as in Sect. f l2.4p is applied, many of 
these pixels are absent. 

Finally, it is interesting to compare monthly mean AOD maps from 
MODIS and TCAM. In Fig. H] the months June and December 2004 are con- 
sidered. Within this figure, TCAM maps account for both moisture and cloud 
mask. In June, while satellite AOD is a factor about 2 larger than model 
AOD, several spatial patterns are recognized in both maps: the eastern Po- 
Valley maximum, the Alps and Apennines minima, the local maximum over 
the Gulf of Genoa, the cloud coverage over western and central Alps. In 
December, AOD values are lower than in June and the departures between 
satellite and model are smaller. However, the number of missing values in 
satellite pictures is much larger than in model computation. This could be 
due to two reasons: i) snow coverage over Alps; ii) fog in the Po Valley. 
While model AOD is delivered in both cases, both meteorological situations 
lead to failure of the MODIS retrieval algorithm. 
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Figure 4: AOD mean maps for June (upper row) and December (lower 

row), 2004 over the QUITSAT domain. Left panels: MODIS product "Cor- 
rected-OpticaLDepth_Land_And_Ocean" from both Terra and Aqua overpasses; Right pan- 
els: TCAM AOD averaged over the times corresponding to the satellite overpasses ac- 
counting for moisture and cloud mask. Notice variable color scales. 
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3.4- Comparisons of timeseries 

Comparisons of spatial maps of AOD from satellite and AOD from the 
model can be useful for checking the ability of the model to reproduce spatial 
patterns. However, satellite AOD itself is affected by sign ificant unce r tainty : 
Minimum relative uncertainty on AOD is as large as 15% iLevy et al.l ( 120071 ) . 
Thus, satellite information alone can not be used as a benchmark of model 
optical performance. The global network of sun-sky photomet ers AERONET 



is often used as a standard reference for AOD measurements Dubovik et al. 



(120001 ). There are 4 AERONET sites within Po- Valley, namely: Ispra, Mod- 
ena, ISDGM-CNR, and Venise. Version 2, level 2.0 AOD values at 870nm 
and 440nm from these sites have been used for estimating AERONET AOD 
at 550 nm. Angstrom power law h a s been assumed to hold in this spectral 
range, as it is done in iHodzic et al.l (l2004f ) 

Comparisons of spatial maps of AOD from satellite and AOD from the 
model can be useful for checking the ability of the model to reproduce spatial 
patterns. However, satellite AOD itself is affected by sign ificant unce r tainty : 
Minimum relative uncertainty on AOD is as large as 15% iLevy et al.l ( l2007l ). 
Thus, satellite information alone can not be used as a benchmark of model 
optical performance. The global network of sun-sky photomet ers AERONET 



is often used as a standard reference for AOD measurements Dubovik et al. 



(120001 ). There are 4 AERONET sites within Po- Valley, namely: Ispra, Mod- 
ena, ISDGM-CNR, and Venise. Version 2, level 2.0 AOD values at 870nm 
and 440nm from these sites have been used for estimating AERONET AOD 
at 550 nm. Angstrom power law h a s been assumed to hold in this spectral 
range, as it is done in iHodzic et al.l (120041 ). 

Comparisons of AOD daily mean values at one of the Po- Valley AERONET 
sites are presented in Fig. [51 Here, daily mean values correspond to an average 
over the times (within 60 min.) of both Terra and Aqua satellite overpasses. 
The Venise site has been selected for its representativeness: it is located at 
the east end of Po- Valley over an offshore oceanographic tower. Instead, both 
Ispra and Modena sites are within a few tenth km from mountains. Finally, 
ISDGM-CNR photometer is located within the city of Venice and is possibly 
affected by local pollution sources. 

In Fig. |5l both whole-2004 monthly-averaged timeseries and zooms on the 
month of June 2004 are shown. Satellite and AERONET information is 
identical in the three rows of the Figure. First of all, one notices that there 
is a fair agreement between MODIS and AERONET AOD, with a maximum 
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Figure 5: AOD timeseries at the QUITSAT cell next to AERONET site "Venise". Left 
panels refer to whole year 2004 while right panels are zooms on the month June 2004. R^d 
lines for AERONET AOD, blue crosses for MODIS AOD, and green circles for TCAM 
AOD. First row: dry; Second row: wet; Third row: wet with cloud mask. 
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monthly value of about 0.4 in Octo ber. For such a high AOD, MODIS 



expected accuracy is better than 30% iLevy et al.l (120071 ). 



Passing to model outputs, in the first row, timeseries of TCAM dry AOD 
are considered. Their values are systematically lower than both satellite and 
AERONET AOD by nearly an order of magnitude. The June-zoom shows 
however that the model follows the same trend of AERONET in the first 
half of the month: in particular, a seemingly exponential build-up of AOD 
is observed from June 6 to June 10. This time interval corresponds to stable 
weather conditions after rains which have interested Po- Valley on the first 
days of the month. 

In the second row of Fig.El TCAM wet AOD is displayed. There is an im- 
pressive model AOD increase with respect to the dry case. In correspondence 
of winter months, the model delivers AOD values which are even higher than 
AERONET ones. Furthermore, minimum values of TCAM AOD are now 
much closer to observations. 

In the third row, the cloud masking mechanism illustrated in Sect. (12.41) 
has been employed. As a consequence, model AOD is nearly halved for 
the first part of year 2004 and the number of June AOD overestimations is 
reduced. The October peak is very well reproduced by both AERONET, 
MODIS, and TCAM. 

4. Conclusions 

An optical module for chemistry transport model TCAM has been de- 
veloped and resulting AOD has been compared to observations over the Po- 
Valley (Italy) during year 2004. 

The optical module takes aerosol concentrations as an input and delivers 
AOD fields as an output. Main inorganics (nitrates, ammonium, sulfates) 
and organic chemical compounds are considered. Within TCAM, a detailed 
description of the aerosol size distribution is achieved through 10 size bins, 
whose position and size can vary during simulations. For the computation of 
AOD, aerosol particles are modeled as homogeneous externally mixed spheres 
and the effect of moisture is accounted for through a parametrization. 

Thanks to the detailed representation of granulometric classes imple- 
mented into TCAM, it is found that the largest contribution to AOD arises 
from inorganic compounds in the sub-micrometric range. This is fully con- 
sistent with comparable literature: Hodzi c et. al. too found th at over 89% of 



AOD arises from the accumulation mode lHodzic et al.l (120041 ): Furthermore, 
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in the WRF-chein modeling exercise, 6 out of 8 size bins whe re centered on 



sub-micron values of the aerosol radius IChapman et al.l (J2009f ) . 

TCAM AOD has been compared to both MODIS AOD operational prod- 
uct and to AERONET AOD. The model is able to capture main spatial and 
temporal features of the observations. However, dry model output is up to 
one order of magnitude smaller than observations. Quantitative agreement is 
significantly improved by the parametrization scheme used for dealing with 
moisture. 

With respect to the issue of quantitive agreement between model and 
observations, one should keep in mind that the final purpose of AOD com- 
putation from TCAM aerosol concentration is to able to assimilate satellite 
AOD into TCAM. Departures between satellite and model AOD represent 
the actual information which can hopefully improve model performance in 
terms of aerosol concentration fields. However, it is important to stress at 
this place that even satellite AOD data are affected by significant uncertain- 
ties. MODIS retrieval algorithm over land implies a relative unc ertainty on 



AOD between 15% and 100%, increasing for lower values of AOD lLevy et al 



( 120071 ). Thus, an assimilation procedure to come has to account for errors of 
both model and satellite AOD in a proper way. To the extent of the compu- 
tation of the optical properties, it is important that the physical features of 
the chemical and transport model are accounted for faithfully. These features 
include both the microphysical aerosol parameters (cp. Tab.[T]), the modeling 
of the extinction coefficient (cp. Eq. ([3])), and the parametrization of moisture 
(cp. Eq. dH])). The best possible implementation of these issues leads to the 
"true" model AOD and corresponding "true" departures between satellite 
and model AOD. In no case, except the case of accidental agreement, such 
true departures can be null, since the model is known to fail in satisfactory 
reproduction of observed PMIO data. Assimilation of such true departures 
into the model will lead to consistent corrections to the model fields. 

The optical module can be used for data assimilation of any satellite AOD 
into TCAM. In particular, the present results suggest that this operation can 
be meaningful already at mesoscale and not only at global and regional scale 
as already tested in the literature. This is in the spirit of QUITSAT project, 
whose aim is to bring together model outputs and observations for a better 
evaluation of air-quality related variables. 
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